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^__< Abstract. We explore the minimum energy configurations of repelling parti- 

f ^ '. cles distributed over n possible locations forming a toric grid. We conjecture 

that the most energy-efficient way to distribute n/2 particles over this space is 
to place them in a checkerboard pattern. Numerical experiments validate this 
conjecture for reasonable choices of the repelling force. In the present paper, 

Uwe prove this conjecture in a large number of special cases — most notably, 
when the sizes of the torus are either two or multiples of four in all dimensions 
C J and the repelling force is a completely monotonic function of the Lee distance 

* between the particles. 

a 

i— i 1. Introduction 

t-H The problem of distributing particles uniformly over a metric space has proved 

an inspiring question for mathematicians, and has also attracted the attention of bi- 
ologists, chemists, and physicists working on crystallography, molecular structures, 
and electrostatics. 



OO 
O 

f^ We consider an interacting particle system with a repelling force between the 

_^J particles. Under influence of the force, the particles arrange themselves in the 

f-) most energy-efficient manner (one that minimises the total amount of energy in the 

£Nj system). The rationale is that repelling particles tend to spread apart until there is 

t-H as much distance between them as possible, which results in some sort of "uniform 

£>. distribution" of the particles. 

• *~j We investigate the case in which the repelling force between any pair of particles 

rS is given by a function / of their distance. This function / can take the form 

& f( x ) — X~ a , with a > 0. This is a generalisation of the inverse square law that 

is well known in physics, where it describes the behaviour of forces like electrical 

charge. It also contains the reduction in power density of an electromagnetic wave 

as it propagates through space, known as path loss. Path loss plays a crucial role in 

the design of wireless networks, for which the exponent a normally ranges from 2 to 

4 (where 2 is for propagation in free space, and 4 is for relatively lossy environments 

[8]). Note that, since / is a decreasing function, we assume that the repelling force 

becomes stronger when the particles move closer to each other. 
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Denote the locations of the N particles by Xi, • ■ ■ , xn- The /-energy experienced 
by particle Xi is defined as ^2i^f(8(xi,Xj)), where 5 denotes the metric. Our 
original motivation stems from the area of wireless networks, in which particles are 
mobile users, and the /-energy of a user is its interference with all other users. 
When experiencing too much interference, a user is unable to transmit data [HE]. 
For a given maximal sustainable interference level T, a natural question is to find 
the maximal packing of users that can be active simultaneously without violating 
the threshold T, i.e., while maintaining /-energy at most T for each user. In the 
case that we will deal with, this question is intimately related to the problem of 
minimising the total energy among all pairs of particles; see Section [2j 

A lot of research in energy minimisation focuses on spheres. In dimension one, 
the uniform distribution on the circle is simply given by the N equally spaced 
points (roots of unity in the complex plane), but in higher dimensions it becomes 
less obvious on what type of patterns the particles arrange themselves, and this 
presents al sorts of mathematical challenges [H [5] . 

In this paper, rather than spheres or other continuous metric spaces, we consider 
a (i-dimensional grid as possible locations for the particles. Further, to avoid bound- 
ary effects, we wrap around the boundaries. This grid has a natural embedding on 
a d-dimensional torus, see Figure [2] which is why we call it the d-dimensional toric 
grid. 

Let Hi be the size of the grid in dimension i, so that we have a total of n = 
n\ ■ ■ ■ rid possible locations for the particles. It seems reasonable to expect that the 
best way, in terms of minimal energy, to distribute n/2 particles over this space 
is to place them in a checkerboard pattern, by which we mean a pattern in which 
exactly one particle is placed at any two neighbouring points. Of course, because 
of the wrapping around, such a pattern exists only when all rn are even. 

For any toric grid, regardless of its dimensions, there are only two such configu- 
rations (see Figure [I]), and the checkerboard pattern is a natural "uniform distribu- 
tion" on a toric grid. It is also a maximal-size independent set in the graph whose 
vertices are the possible locations and whose edges connect locations differing by 
±1 in one coordinate and not at all in the remaining coordinates. The Lee distance 
between two locations is defined as the length of the shortest path, in terms of 
edges, between the corresponding vertices on this graph. The following conjecture 
formalises, to some extent, our intuition that checkerboard patterns are optimal. 

Conjecture 1.1 (Checkerboard Conjecture). For any toric grid with m, . . . ,rid all 

even, and for any reasonable function f of the Lee distance, the arrangements of 
N = n/2 particles that minimise the maximal f -energy experienced by any of the 
particles are the two checkerboard patterns, and only these. 

Of course, the validity of this conjecture depends on what reasonable means. It 
will certainly imply decreasing, but more conditions will be necessary. 

For the example of wireless networks the Checkerboard Conjecture has the fol- 
lowing implication. When the threshold T is increased slowly, an increasing number 
of users is allowed to be active simultaneously, and there will be a critical value 
T = T* at which this number becomes precisely n/2. At that critical threshold, we 
then know that n/2 active users can only be arranged in a checkerboard pattern — a 
phenomenon that we have first observed in simulations of wireless networks. For a 
large toric grid it is of course infeasible to prove the Checkerboard Conjecture by 
exhaustive search through all (^) possible configurations. 
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Figure 1. A 6 x 6-grid, the black and white points denote the 
checkerboard configurations. 




Figure 2. An embedding of a 6 x 6-grid on a 2-dimensional torus. 

In this paper we prove the Checkerboard Conjecture in a large number of special 
cases. We first derive a continuous relaxation of our energy minimisation problem 
in Section [2j Then, as a first (and known) special case, we prove the Checkerboard 
Conjecture for the one-dimensional toric grid and f(x) — 1/x in Section k5j In 
Section [3] we prove the Checkerboard Conjecture for the toric grid that has size two 
in any dimension and functions / whose fc-th forward difference has sign (— 1) . 
Section ^ then gives our strongest result (which however does not immediately 
imply any of the previous two results): a proof for the Checkerboard Conjecture 
in the case where each n, is either two or a multiple of four, now assuming that 
/ is completely monotonic; see Section [5] for this notion. Finally, in Section [6] we 
discuss some problems open for future research. 



2. Framework 

In this section we formulate a general version of our energy minimisation prob- 
lem, which belongs to the realm of discrete optimisation. We then derive a continu- 
ous relaxation of that problem. If all optimal solutions of the relaxation happen to 
be feasible for the discrete problem, then they are also the optimal solutions for the 
discrete problem. This is our strategy for proving special cases of the Checkerboard 
Conjecture in the following sections. 

Let d be a positive integer, the dimension of our toric grid. Let m, ■ ■ ■ , n<j 
be positive integers, and set G := Ili=i(^/ n i)7 a finite Abelian group. Let A := 
{(Si 9) I 9 G G} b e the diagonal in G x G. Let u be a function (GxG)\A->t that 
is symmetric (u(g, h) = u(h, g) for all g,h £ G) and G-invariant (u(k + g,k + h) = 
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u(g, h) for all g, h, k G G). Later we will add the restriction on u that it is a suitable 
(and in particular decreasing) function of some G-invariant distance on G, so that 
we may think of u(g, h) as a repelling force between particles located at g and h. 
But for deriving our relaxation this restriction is not needed. 

Next let p be a natural number less than or equal to \G\. Let S be a subset of G 
of cardinality p; this is the set of locations of our particles. Given h € S we define 

E h (S):= J2 M ^5)< 
ges\{h} 

called the energy of S experienced by h, and 

E max (S) := ma,xE h (S), 

called the maximal energy of S. The Checkerboard Conjecture concerns the follow- 
ing optimisation problem with p = |G|/2. 

Problem 2.1 (Maximal energy minimisation). Minimise the maximal energy E max (S) 
over all subsets S C G of cardinality p. 

Further we define 

E tot (S):=J2 E h(S), 
hes 
called the total energy of S. We can now formulate a second optimisation problem. 

Problem 2.2 (Total energy minimisation). Minimise the total energy E to t{S) over 
all subsets S C G of cardinality p. 

An important relation between these two optimisation problems is the follow- 
ing. If a set S* minimises the total energy and happens to have the property 
that Eh(S*) = E g (S*) for all h,g <G S* , then S* also minimises the maximal 
energy — indeed, any S with -E ma x(<S') < E max (S*) would necessarily have E tot (S) < 
E tot (S*). The condition that Eh(S*) does not depend on h e S* is satisfied, in 
particular, when S* is a coset of a subgroup G' of G. Indeed, the set S* \ {h} is then 
equal to G' for all h e S* , and hence Eh{S) — J2 g eG' u (^' ^ + 9) = J2 g eG' M (0> #)' 
independently of h € S* , by the G-invariance of u. Checkerboard patterns S are ex- 
amples of coscts. In view of this relation between maximal energy minimisation and 
total energy minimisation, most of our paper focuses on total energy minimisation. 

We proceed to derive a continuous relaxation, which uses the space V = K 
of real-valued functions on G. We represent a subset S C G by its characteristic 
vector xs € V defined by 

. . ) 1 if g € S, and 
%s{g) = < n .„ . „ 
[0 iig£S. 

The space V has two important additional structures. First, it is equipped with a 
G-action defined by (gx)(h) — x(h — g) for g,h G G and x e V; here the minus 
sign is customary but not strictly necessary. Second, it has a natural inner product 
(x\y) := ^2„ e Qx(g)y(g) that is G-invariant in the sense that (hx\hy) = (x\y) for 
all h€G. 

Now we define the key quantity in this paper: the energy kernel 

A -= J2 "(5,0)3, 
gec,g^o 
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which lives in the group algebra RG and acts by a linear map on V: 
(Ax)(h) = J2 u(9,0)x(h-g). 

Alternatively, one may think of A as the symmetric G x G-matrix whose (g, /i)-entry 
equals u{g, h), and perhaps this conceptually simplifies the following computations. 
For any subset S of G we have 

(x s \Ax s ) = J2 x s(h)(Axs)(h) = J2(Ax s )(h) =J2 J2 u(g,0)x s (h- g) 

h£G heS h£.Sg£G,g^0 

= E E u(h-g,0)x s (g) = Y^ E <h~9,0) 

h£Sg£G,g^h heS geS,g^h 

= E E u ( h >9) = E tot (x s ). 

heS geS,g^h 

Thus we have identified the total energy of a set S as the inner product of xs 
with its image under the energy kernel. Moreover, if S has cardinality p, then xs 
has the further properties (xs\xs) = (xs|l) = p, where 1 is the all-one vector in 
V. This motivates the following optimisation problem. 

Problem 2.3 (Fractional energy minimisation). Minimise the "fractional total 
energy" (x\Ax) over all x £ 7 satisfying the constraints (x\x) — (x\l) = p. 

Note that if we add the additional constraint that x be a 0/1-vector, then this 
problem is just total energy minimisation. Thus a feasible solution of the relaxation 
can be thought of as a "fractional" solution to the original problem, though this is 
slight abuse of language for two reasons. First, the optimal solution to the problem 
just stated may very well have irrational coordinates, so that it is not a "fraction" . 
Second, it may have coordinates that are not in the interval [0, 1], so that it cannot 
be "rounded" to an actual solution of the original problem. We could have added 
the additional constraint that all coordinates lie in that interval. This would have 
led to a more general quadratically constrained quadratic program. However, as 
we describe next, the formulation just given allows for a simple solution once the 
eigenvalues of A are known. But first we discuss an example to make all notions 
more concrete. 

Example 2.4. Let d = 2 and m = n 2 = 4, so that G = (Z/4) x (Z/4). We depict a 
vector x € V as a 4 x 4-matrix of real numbers, with rows and columns labelled by 
the elements 0, 1, 2, 3 of Z/4 and entry x(k, I) at position (fc, I). Here is an example 
of the action of G on V ( "shift 1 down and 2 to the right" ) : 



(1,2) 



12 3 4 

5 6 7 8 

9 10 11 12 

13 14 15 16 



15 16 13 14 

3 4 12 

7 8 5 6 

11 12 9 10 



In this case the Lee distance 5 on G is given by 

5((j,k),(l,m)) = \l-j\ + \m-k\ 
where \a\ is the smallest non-negative integer in (a + 4Z) U (—a + 4Z). Define 



u{g,h) 



S(g,h) 



, so that particles situated at g and h repel each other with a 
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force inverse proportional to the Lee distance between g and h. Take p — 4 and 
S = {(0, 0), (0, 1), (0, 2), (0, 3)}. Then the energy kernel A satisfies, for instance, 



Axs = A 



1 


1 


1 


1" 







































+ 1 



2 ' 3 

i + 1 



+ 1 



i ! 



1 + 5 + 1 

1,1,1 



and we find that (x|Ax) equals 4(0 + 1 + 1/2 + 1) = 10, which is indeed the total 
energy of S. Note that one could represent A by a symmetric 16 x 16-matrix. 

The geometric intuition behind what follows is that, to minimise (x|^4x) over 
the sphere defined by the constraints (x|x) = (x|l) = p, we have to maximise the 
component of x in the eigenspace of A corresponding to the smallest eigenvalue of 
A. Indeed, since A "is" a symmetric matrix, or, equivalently, a self-adjoint linear 
map with respect to (.|.), we know beforehand that all of its eigenvalues on V are 
real. Nevertheless, to determine those eigenvalues it is convenient first to complexify 
V to Vc := C <£) V = C G , because this allows us to simultaneously diagonalise all 
group elements g in their action on V. As a consequence, their linear combination 
A is then also diagonalised. 

We use basic terminology from the representation theory of finite Abelian groups, 
for which we refer to [TTJ Part 1]. Since G is Abelian, Vc splits as a direct sum of 
simultaneous eigenspaces of all g £ G, that is, we have 

^c = V x , 

xec v 

where \ runs over the set G v of group homomorphisms (G, +) — > (C, •) (called the 
character group of G) , and where V x is the subspace of V defined as 

V x := {v € V c | gv = x(g)v for all geG}. 

Distinct V x are orthogonal with respect to the complexification of (.|.) to a Hermit- 
ian inner product, and as Vc is isomorphic to the regular representation of G, each 
V x is one-dimensional. More explicitly, being a group homomorphism, \ € G v is 
determined by its values Ci, • • • ,C,d on the generators (1, 0, . . . , 0), . . . , (0, . . . , 0, 1) of 
G. Each ^ is an n^-th root of unity, and the space V x is spanned by the element of 
Vc whose value at (g\, . . . , gd) equals (i 91 ■ ■ ■ (^ 9d . Note that this function is just 
the multiplicative inverse x _1 of x, which equals the complex conjugate \ since the 
C,i lie on the unit circle. 

Since A is a linear combination of the group elements, it acts by scalar multipli- 
cation on each V x , and the scalar equals 

A(x):=5>(<7,0)x(<?). 

9#0 
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Note that this is the (discrete) Fourier transform of the function G — > C, 
u(g, 0). As noted before, this is a real number, since 

W)=Y1 <9, o)xig) = E <9, o)x(-g) = A( x ), 

9#0 ff #0 



where the second equality follows from the fact that x(#) = x(d) = x(~ #)j an d 
where we have used the symmetry and G-invariance of u in the last step. This 
computation also shows that X(x) equals A(x). Thus the real vector space V splits 
as a direct sum 

V = © ^ { x,x} 
{x,x> 
over unordered pairs of a character and its conjugate (which may coincide), where 
V{x-x} ^ s * ne rea, l s P ace °f real-valued vectors in the complex space V x + V%. Here 
the direct sum is in fact orthogonal, each summand is one-dimensional or two- 
dimensional according to x being real or non-real, and each summand has corre- 
sponding A-eigenvalue A(x). As a special case, the trivial character 1 sending every 
element of G to 1 is real, and the vector space Vm is spanned by the all-one vector 
1. Now we can solve fractional energy minimisation as follows. 

Proposition 2.5. Let \ m i n be the smallest value of \(x) «s X ranges over the non- 
trivial characters. Then the set of optimal solutions of fractional energy minimisa- 
tion consists of all vectors of the form rgr 1 + y where y belongs to the eigenspace of 
A inV with eigenvalue X m in, is perpendicular to 1, and satisfies (y\y) = p — p 2 /\G\. 

Note that the condition that y be perpendicular to 1 is automatic if A m j n < A(l), 
because eigenvectors corresponding to distinct eigenvalues are perpendicular. In 
this case the set of optimal solutions forms a sphere whose dimension equals the 
multiplicity of A m i n minus one, embedded in the affine hyperplane where (x|l) = p. 
If A m i„ happens to equal A(l), then the dimension of that sphere of optimal solutions 
equals the multiplicity of A m i n minus two. 

Proof. Consider any feasible solution x. The condition (jc|1) = p means that x can 
be written as t^tI plus y with y in the direct sum of all Vr^v} with y non-trivial. 
The condition that (x\x) = p translates in the above condition on the norm of y. 
Let y{ x ,x) be the component of y in Vi x ^\. Then, by orthogonality, 

11 ' ' ' ' {X:X> non-trivial 

For this expression to be minimal among all y with squared norm p — p 2 /\G\ it is 
necessary and sufficient that all components yi x ,x} f° r which A(x) is not equal to 
Amin are zero. This proves the proposition. □ 



(x\Ax) = (-^1 + y\A(£-l + y)) = A(l)£- + ]T A(Y)|| y{x ^ } || 2 . 



Example 2.6 (Continuation of Example |2.4[ ). In this case a character x G G v is 
determined by its values Ci := x(l, 0) and C2 := x(0j 1): an d these complex numbers 
must satisfy Ci = C2 = 1- Conversely, any such pair determines a character, which 
maps the pair (51,52) G (<^/4) 2 to CiXf 2 - The eigenvalue A(l,i) of A equals (reading 
elements row-wise) 

1 . . i 1 i , ,1 i 1 is, i 1 * x 13 

+ i i) + (l+ ) + (- + + 1 + ) = — • 

v 2 ' v 2 3 2' v 2 3 4 3' v 2 3 T 12 

The eigenvalues of the energy kernel A are as follows. 
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12-A(Ci,C 2 ) 
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This shows that A(l, 1) is the maximal eigenvalue, with a one-dimensional eigenspace 
spanned by 1; and A(— 1,— 1) is the minimal eigenvalue, with a one-dimensional 
eigenspace spanned by 



Z = ((91,92) -► (-1)* . (-1)») = 



Hence fractional energy minimisation has two optimal solutions (forming a zero- 
dimensional sphere). For p = 4 these are |l± 4 #• Both of these have coordinates 
in [0, 1], but neither of these is integral or even rational. Of course, a brute-force 
computer search over all configurations of p = 4 particles on G is possible in this 
case. Such a search reveals that, up to translations, the following three character- 
istic vectors are the ones minimising the total energy: 
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Note that the first two are each other's images under swapping rows and columns, 
and that they represent the cyclic subgroups of G generated by (2, 1) and (1,2), 
respectively. The last one is not a coset of a subgroup, but still has the property 
that the energies experienced by all particles are the same. This means that these 
configurations also minimise the maximal energy. We conclude that for p = 4 the 
optimal solution to the fractional energy minimisation problem has little bearing on 
the original discrete optimisation problem, and it would be very interesting to find 
a strengthening of our relaxation that does give the optimal solutions for p = 4. 

On the other hand, consider now p = 8 = \G\/2. Then the optimal solutions for 
fractional energy minimisation are jlljZ, and these are exactly the characteristic 
vectors of the two checkerboard patterns. Since fractional energy minimisation is a 
relaxation of total energy minimisation, this proves that the checkerboard patterns 
also minimise total energy -Etot- And since they are cosets of a subgroup, the 
argument given after the introduction of E to t shows that the checkerboard patterns 
are also the (unique) minimiscrs of E max . We have thus proved our first instance 
of the Checkerboard Conjecture. In the following three sections we follow exactly 
the same strategy: we determine the minimal eigenvalue of the energy kernel to 
determine the optimal solutions of fractional energy minimisation, and when these 
are the checkerboard patterns, then we are done. This works in quite a number of 
cases — in particular, when each rii is either two or a multiple of four, see SectionlS] — 
but not in general, see Section [6j 

We conclude this section by relating our techniques and results to existing liter- 
ature on energy minimisation. The most relevant reference is [1], which is primar- 
ily concerned with energy minimisation of repelling particles on high-dimensional 
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spheres. There are close parallels between our approach and the techniques there. 
In particular, both papers make essential use of the representation theory of the 
natural symmetry group of the problem. In our case this is the group G itself, 
while in the spherical case it is the orthogonal group. The fact that G is Abelian 
(and hence diagonalisable) makes the representation theory much easier than that 
required for the spherical problem. On the other hand, there are also important 
distinctions. First of all, [I] makes essential use of the fact that the symmetry 
group acts distance-transitively, while our group G acts only simply transitively 
on itself. Second, while [4 deals with continuous optimisation problems, our op- 
timisation problem is inherently of a discrete nature. It could be strengthened 
into a continuous problem by letting " 1 ' 2 '" d particles move on the continuous torus 
R d /(n]Z x • • • x n<jP). If the energy-minimising configurations for that problem are 
(unique up to translation and) equal to our checkerboard configurations, then that 
implies our results (provided that one uses the 1-norm, which is the continuous 
analogue of the Lee distance). It would be interesting to know if the techniques 
from [4] can be used to prove such a stronger statement. Note that the pre-image 
in M. d of the particles' positions form a periodic set. The paper [3] does discuss op- 
timisation problems for such sets, but in the context where the lattice of periodicity 
is allowed to vary, as well. In our case, the lattice is fixed to n\L x • • • x n<iL. 

Our minimisation problem is also reminiscent of certain instances of the quadratic 
assignment problem that are studied in [3]. In that paper it is proved that the class 
of those instances is NP-complete. In view of this, minimality of checkerboard 
configurations is an interesting result. On the other hand, in [3] the corresponding 
maximisation problem is proved to be solvable in polynomial time. In our setting, 
this would correspond to taking for u some increasing function of the distance, 
which has an interpretation in terms of attracting particles rather than repelling 
ones. This problem, in turn, resembles the problem of how to build an optimal 
city [2]. Our attention, however, is focused entirely on the setting where u is some 
decreasing function of the Lee distance. 

3. The one-dimensional case 

The first special case that we consider is the one-dimensional toric grid. That 
is, we consider G = (Z/n), with n a positive integer, and define a graph on G 
by connecting elements that differ by ±1. One can think of this as n equally 
spaced points on a circle. Further we define the function u: GxG\A-> M + 
as u(g,h) :— f(6(g,h)), where S(g,h) denotes the Lee distance between g and h 
and / : M.+ —> R + is a decreasing strictly convex function. In [Fj the continuous 
problem of letting particles move on a circle is considered and it is shown that total 
energy is minimised if and only if the particles are equally spaced. This validates 
the Checkerboard Conjecture for this special case. 

The proof in [5] is based on a direct comparison between the total energy in the 
configuration with equally spaced particles (i.e. the checkerboard configuration and 
its translations) and the total energy in any other configuration. To get a better 
feeling for the technique of Section^ we now assume f(x) = x~ x and prove that 
fractional energy is uniquely minimised as well by the checkerboard vectors. That 
is, we will prove the following result. 

Theorem 3.1. Assume that n is a multiple of 2. Then the fractional energy (x\Ax) 
over all x € V = M. G with (x\x) = (x\l) = n/2 is minimised by the characteristic 
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vectors Xs mea and xs odd of the checkerboard configurations 
Seven/odd ■= {g € G \ g even/odd}, 

and only by these. As a consequence, both the total energy and the maximal energy 
among all subsets S C G with \S\ = n/2 are uniquely minimised by S even and S dd- 



By Proposition [^5] we know that set of optimal solutions of the fractional energy 
minimisation problem is determined by the smallest eigenvalue of the energy kernel 
A. As discussed in Section [2] the eigenvalues are functions on the dual group 
G v = C n , with C n the group of complex n-th roots of unity. And the fractional 
energy is uniquely minimised by the checkerboard vectors if the eigenvalue has 
a unique minimum at —1, which lives in G v as n is even. Thus the following 
proposition implies Theorem |3.1| 



Proposition 3.2. The function X : G v -¥ R, 

a(x) = £«%.())- V 

has a unique minimum at x — ~ 1 ■ 
Proof. As n is even we can write 

a/2-1 

1 _;. . i- . - 

x 



A (x) = E \(x- k + x k ) + 2 r n - 
fc=i 



Now note that the n elements of the group C n can be represented by e 27ri i/ n for 
j € J = {0, . . . , n — 1}. With a slight abuse of notation we can thus write 

" /2_1 1 f 2njk, 2 



m= E ^cos(^) + - ( -ir 

' ' h mm 



k=l 



So, we need to prove that argminA(j) = ^. 
For this first define k(x) : [0,n) — >• K by 

n/2-l 



L z — ' k n n 

fc=i 



Now note that all eigenvalues are real numbers, so n(j) = A(j) for all j e J . Hence, 
argmin k{x) — § implies argminA(j) = ^. 

To prove this let us look at the derivative of k(x). Differentiating k(x) we get 



n/2-1 
-47T v-^ . ,2TTXK S Z7T 



k'(x) = -^ Y^ sH^r 1 ) 



smixn) 



fc=i 
To simplify this expression note that 



n 



cos(f) — cos((m + \)x) 



( sin(fcx 



^- ' 2sin(f) 



as can be proved using the Prosthaphaeresis formulas. 
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We thus find 

« (*) = " , t— r — sm(7rx) 

v ; n 2 sin ^) n v ; 

-4ttcos(^)- (cos(7Tx)cos(^) + sin(7ra;)sin(^)) 27T . . , 

■ sm(7rx) 



n 2sin(^) n 

= — COS 7Tx) - 1 , " . 

n sm — 

^ n / 

Hence, k'(x) < for x e (0, | ) and «/(x) > for x e (§, ra). Further, by l'Hopital's 
rule, lim x _y K'(x) = 0. Also, k'{x) < for x £ (f — 1» f ) and «'(x) > for 
x £ (tJ, tt + 1). Thus k(x) has a unique minimum at x* = ~, which completes the 
proof. D 

4. The Hamming cube 

In the previous section we considered the one-dimensional toric grid with an 
arbitrary number of points. In this section we consider the d-dimensional toric 
grid, but now we will assume that the grid has size two in all dimensions. That 
is, we consider G := Ili=i(^/2) — {0, l} d - Further we define a graph on G by 
connecting two elements if they differ by (0, . . . , 0, 1, 0, . . . , 0) for some position of 
the 1. One can think of this graph as the ci-dimensional unit hypercube; the corners 
of this hypercube are the points of the toric grid. In fact, one should construct a 
d-torus from the d-dimensional hypercube, but this only connects corners that are 
connected already, as our grid has size two in every dimension. 

Let / : Z + — > R + be such that its m-th forward difference has sign (—1)™, i.e., 
we assume that 

(1) (-l) m A™[/](x)>0 

for all x 6 Z + . Here A denotes the forward difference, A[/](x) = f(x + 1) — f(x). 
We now define the function u : G x G \ A — > R + as u(g, h) := f(S(g, h)), where 
8{g, h) denotes the Lee distance between g and h. We will validate the Checkerboard 
Conjecture for this case. That is, we will prove the following theorem. 

Theorem 4.1. The fractional energy (x\Ax) over all x £ V — M. G with (x\x) = 
(x\l) = \G\/2 — 2 d ~ l is minimised by the characteristic vectors xs emn and xs odd of 
the checkerboard configurations 

Seven/odd ■= {(ffi, ■ ■ ■ , 9d) & G | ^2 g l even/odd}, 

i 

and only by these. As a consequence, both the total energy and the maximal energy 
among all subsets S C G with \S\ = 2 d ~ 1 are uniquely minimised by S even and S dd- 



To prove Theorem |4.1| we again capitalise on Proposition |2.5| and determine 
the smallest eigenvalue of the energy kernel A. As the toric grid has size two in all 
dimensions we know by the discussion in Section[2]that the eigenvalues are functions 
of the dual group G v = { — 1, l} d . We further know that the fractional energy is 
minimised by the checkerboard vectors if the eigenvalue has a unique minimum at 
(— 1, . . . , — 1) £ G v . This we prove in the following proposition, which thus implies 
Theorem l4Tl 
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Proposition 4.2. The function X : G v -» R, 

d 

A(x)=E/(^,0))nxf i 

/las a unique minimum at x — {— lj ■ • ■ , — 1)- 

Proof. Because of symmetry it suffices to prove that if x, X G {~ 1; l} rf are characters 
that differ only in the first position, where xi — 1 ; Xi = — 1) then A(%) > A(x). Wc 
therefore compute 

A(X) - A(X) = E /(%> ))(! - (-!) 9l )xr • • • Xf 

= 2 y, /(%,o) + f) X f •••xr 

S=(S2,---,9d) 

Writing g for the number of ones among the X2> • ■ • > Xd and using dummy I for the 
number of indices j with Xj = 9j = 1 an d dummy k for the number of indices j 
with Xj = — f and g rj = 1 the above expression equals 

»£(!)sT/<<+*+i>C";"*w 

= 2^(^(-l) d - 1 -«(A d - 1 -V)(« + l) 

z=o ^ ' 
>0; 

where the last inequality follows from the conditions on / and the last equality 
follows from the well-known formula 



k=0 



A m [/](^) = E(r)(- 1 ) m+fc/ ( a; + fc ) ; 



see for instance [TD], page 203, equation la]. This proves the proposition. □ 



5. Multiples of four 

In this section we prove our strongest result about the Checkerboard Conjecture, 
concerning {n\ ■ ■ ■ rid)/ 2 particles on the Abelian group G = (Z/mZ) x • • • x (Z/rirfZ) 
when each rn is either two or a multiple of four. Define a graph on G by connecting 
two elements if they differ by ±(0, . . . , 0, 1, 0, . . . , 0) for some position of the 1. The 
Lee distance is the shortest-path distance between g, h G G in that graph, and 
denoted by 5(g, h). The energy will be measured with a function / : R + — > R + that 
we assume to be (strictly) completely monotonic, i.e., C°° and with (— l) k f( k > (x) > 
for all x € R+ and for all k. 

Remark 5.1. Typical examples of completely monotonic / are f{x) = x~ a for a a 
positive real number. These functions appear prominently in the paper [4] , as well. 
Moreover, as mentioned in the introduction, these function generalise the inverse 
square law in physics, and describe the path loss for electromagnetic waves. 
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Remark 5.2. Note that if g : M+ — > R + is completely monotonic, then so is — Ag, 
because 

nx+l 

(-l) k (-Ag)( k \x) = {-l) k {g (k \x) - g (k) {x + 1)) = (-l) fc+1 / g {k+1) {t)dt > 

J X 

for all x > 0. Iterating, we find that (—A) k g is a positive function, so that the 



restriction of g to Z + satisfies the condition on /. In particular, Theorem 4.1 implies 
the special case of the upcoming Theorem |5 . 3| where all rn are equal to 2. 



Define the function h:GxG\A4 W. + as u(g, h) := f(S(g, h)). The maximal 
energy £ max , the total energy -Etot, the energy kernel A, and the fractional energy 
are defined using this function u. 

Theorem 5.3. Assume that each ofn\, ...,Ud is either 2 or a multiple of A. Then 
the fractional energy (x\Ax) over all x e V = M. G with (x\x) = (x\l) = \G\/2 is 
minimised by the characteristic vectors Xs emn and xs odd of the checkerboard config- 
urations 

Seven/odd •= {(ji, ■ ■ ■ , jd) \ ^Ji even/odd}, 

i 

and only by these. As a consequence, both the total energy and the maximal energy 
among all subsets S C G with \S\ — \G\/2 are uniquely minimised by S even and 

Sodd- 



By Proposition 2.5 it suffices to find the smallest eigenvalue of the energy kernel 
A. Recall that the eigenvalues are functions on the dual group G v , which in our 
present setting we may identify with G v = C ni x ■ • • x C ni , where C a is the group 
of complex a-th roots of unity. Since all n, are even, G v contains the real character 
— 1 = (— 1, . . . , — 1), with character space V{_i} spanned by the vector z defined by 
z(ji, . . . , jd) = ( — iy i +---+3t. The checkerboard configurations ha ve c haracteristic 
vectors |l ± \z. Thus the following proposition implies Theorem 5.3 

Proposition 5.4. The function G v — > R sending \ t° 

has a unique minimum at (—1, . . . , — 1). 

Proof. The completely monotonic function / can be written as 

Jo 
where a is non-decreasing and the integral converges for < x < oo (for this 
theorem of Bernstein see [12 Chapter IV, Theorem 12b]). Thus, after interchanging 
integral and (finite) sum the function that we seek to minimise is seen to equal 



X^ 



J™ ^e-^ X (g)) da(t). 



For this function to have a unique minimum at x = (— 1, . . . , — 1) it suffices that 
for each fixed t > the function 

G V ^M, X ^ 5>-^' 0)t x(3) 

9^0 
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has a unique minimum at x = (— 1, . . . , — 1). At this point, there is no harm in 
including g = in the sum, which adds a constant term 1 independently of \. 
Writing a :— e* > 1, we want to minimise 

£«- 5(9 < 0) x(s) 

9 

over all characters \. Write g = (.gi, . . . ,gd)] then the Lee distance S(g,0) equals 
S(gi, 0) + • • • + 5(gd, 0), where we use the notation d also for the Lee distance in the 
individual Z/n,Z. Hence we have 

9 (fflv.Sd) 

which factorises as the product of the factors 

for i = l,...,d. Thus we are done if we can show that each of these factors is 
positive for all \i, and minimal at Xi = ~~ 1- 

What remains is an easy calculation. Simplify notation by fixing i and writing 
n := rii. As a warm-up, consider the case where n = 2. Then the trivial character 
1 of Z/2Z is mapped to 1 + a -1 , while the character —1 of Z/2Z is mapped to 
1 — a -1 . The latter value is smaller than the former, as desired. Next assume that 
n is divisible by 4, and consider the function A : C n — > R sending an n-th root of 
unity ( to 

J2 a-^>°>C s 



n/2-l 

i+ y a - 9 (c 9 +c 9 )+a^ n/2 c /2 



9=1 



n./2-l 



9=0 



i -n/2/n/2 i „-n/2/-n/2 

= t 4 + ; -I7-I - u - «-" /2 c- /2 ). 

1 — a •% 1 — a 1 ^ 1 

Now C™' 2 equals ±1. If it equals 1, then this expression reduces to 

(1 - '-""> (t^c + t^t " = » - °-" 2 » (iii^w) • 

Since a > 1, this is a positive number for any £, and its smallest value is attained 
for the n/2-th root of unity £o furthest away from 1 in the complex plane, which 
equals —1 since n is a multiple of 4. Next, if (J 1 / 2 equals — 1, then the expression 
above reduces to 

^"Mt^c + t^t^- 1 ) =a + «-' 2 )(p i ^ P ). 

which is again positive for all C, with <^ n ' 2 = — 1 and minimal for the d furthest 
away from 1 in the complex plane (which equals — e ±27 ™/™). 
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Figure 3. Eigenvalues of the energy kernel (interpolated to make 
an easily interpretable surface; and truncated near the corners). 



It remains to compare A(£o) and A(Ci)- For this, note that the factor (1 — a~"/ 2 ) 
in the expression for A(£o) is smaller than the factor (1 + a~ n / 2 ) in the expression 
for A(Ci), while moreover £o is further away from 1 than £i. Hence the minimum is 
attained at £ = £q = — 1, as desired. This concludes the proof of the proposition, 
and hence the proof of Theorem |5.3| □ 



Remark 5.5. The problem in extending this argument to the case of general 
even rii is that in the proof of the proposition, when n/2 is odd and at least 3, 



A(Ci) = A(— 1) can very well be larger than A((o)- In other words, Proposition 5.4 
does not extend to that situation for general completely monotonic functions / 
Extensive numerical evidence with d = 2 and ni,ri2 up to 100 suggests that the 
proposition does extend to the case where f(x) = 1/x (and probably even for 



f{x) — x~ a with a > 0), in which case also Theorem 5.3 extends. 

Of course, the Checkerboard Conjecture may still very well hold for general 
completely monotonic / and general even n^, even when Proposition |5.4| does not. 
To prove this, one would like to have a combinatorial/ algebraic argument why the 
eigenspaces corresponding to characters x having some entries equal to the relevant 
— exp(27ri/n) with n/2 odd are no good for finding 0/1-vectors. At this moment 
we have no such argument. 

Remark 5.6. Note that the table of eigenvalues of A contains \G\ entries, and 
evaluating each individually involves a summation over \G\, so that a straight- 
forward implementation would require C(|G| 2 ) (floating-point) operations. How- 
ever, the table of eigenvalues is nothing but the Fourier transform of the table 
of values f(5(g, 0)), so that algorithms for (d-dimensional) Fast Fourier Trans- 
form may be used to compute the table of eigenvalues more efficiently. We have 
used Mathematical built-in function Fourier in our numerical experiments. For 
instance, Figure [3] shows the eigenvalues of the energy kernel for d — 2 and 
n\ = n 2 = 10 and f{x) = x~ ' 3 as a function of (j,k), standing for the char- 
acter (eT^Tr). The "rough edges" are related to the even/odd effect in the 
proof above, but they smooth out in the middle, so that the minimum is attained 
at j = k = 5, which corresponds to the character (—1,-1). 
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6. Open problems 

In this paper we have proved the Checkerboard Conjecture for toric grids and 
energy functions / for which the unique minimal eigenvalue of the energy kernel A 
is uniquely attained at the character (— 1, . . . , —1) in the dual group. Indeed, the 
minimal eigenvalue of A always governs the solution to fractional energy minimi- 



sation (Proposition 2.5), and if it is uniquely attained at (— 1, . . . , —1), then the 
minimisers of fractional energy minimisation are 0/1-vectors and hence minimisers 
of the original problem. This motivates the following open problems. 



6.1. Generalising Theorem 5.3, One would like to extend the results of Sec- 
tion [5] to the case where some of the rii are larger than 2 but equal to 2 modulo 4. 
Computer experiments with two-dimensional grids up to 100 x 100 suggest that for 
f(x) = 1/x (and perhaps for general inverse power laws) the minimal eigenvalue of 
the energy kernel is still uniquely attained at (— 1, . . . , — 1). However, a proof of this 
fact would require more sophisticated analytical methods than the factoring into 
one-dimensional instances that we used there. Alternatively, one might try and find 
a combinatorial argument why the minimisers of total energy minimisation are nec- 
essarily perpendicular to the eigenspaces of A corresponding to the "undesirable" 
characters, thus sharpening the relaxation. 

6.2. Euclidean or other metrics. The question remains for which metrics other 
than the Lee distance the Checkerboard Conjecture will hold. For the Euclidean 



metric the analogue of Proposition 5.4 does not hold for general completely mono- 
tonic functions /. The first steps of its proof go through when one replaces the Lee 
distance with the square of the Euclidean distance. This holds in particular for the 
all-important factorisation step. One is then lead to minimise 

J2 a-^ 0) V 



over the n-th roots £ of unity, where S stands for the Euclidean metric on the one- 
dimensional toric grid, and where a is a real number greater than 1. However, for 
a close to 1 the minimum is not in general attained at —1; see Figure [4] On the 
other hand, extensive numerical evidence for 2-dimensional grids up to 100 x 100 
suggests that the analogue of Proposition |5 ,4| does hold when the energy is given by 
an inverse power law of the Euclidean distance (and if these computations were done 
in suitable interval arithmetic, they would prove the Checkerboard Conjecture for 
those grids). Perhaps one can prove such a statement by bounding the contribution, 
in the Bernstein expansion of these inverse power laws, of the bit where a is small. 
The Checkerboard Conjecture is not true for every metric however. For example, 
for the Chebyshev (or Loo) metric one can compute that the maximal energy of 
the patterns in Figures 5(a) and |5(b)| is lower than the maximal energy of the 



checkerboard pattern in a 6 x 6-grid. It is unknown what properties a metric 
should have in order for the Checkerboard Conjecture to be true. 

6.3. Nearly optimal configurations. It is also interesting to consider nearly 
optimal configurations of n/2 particles in toric grids. When considering total energy 
E to t , these configurations are probably given by small perturbations of checkerboard 
patterns. However, nearly optimal configurations for maximal energy E max seem to 
exhibit beautiful, regular patterns. For example, for a 6 x 6-grid we have computed 
the maximal energy of all possible ( X g) configurations for the function f(x) = A 
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Figure 4. Values (vertical axis) of ^ & 

exp(^p), k = 1, . . . , 7 (horizontal axis); k = is excluded to make 
it clear where the minimum is attained: at k = 2, 6 rather than 
k = A. 
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(a) Second best for Lee and Euclidean metric, (b) Third best for Lee and Euclidean metric. 
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(c) Fourth best for Lee metric. 
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(d) Fourth best for Euclidean metric. 



Figure 5. The black nodes denote the nearly optimal configura- 
tions for £ max in a 6 x 6-grid for f(x) = -. 



and for both the Lee metric and the Euclidean metric. This computation shows 
that the Checkerboard conjecture holds (a fact that one can prove more efficiently 



using the discrete Fourier transform as in Remark 5.6 ), but also that the second and 



third best configurations agree for the Lee and Euclidean metrics (see Figures 5(a) 



and 5(b) ). The fourth best configurations are different for the two metrics though, 
see Figures 5(c) and |5(d)| It is a widely open question how to prove that certain 
configurations arc second or third best. 
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